The Energy Momentum Tensor of Chiral Quark Models at Low Energies 



E. Megi'asQ E. Ruiz ArriolaQ and L. L. Salcedc0 

Departamento de Fisica Moderna, Universidad de Granada, E-18071 Granada, Spain 

(Dated: February 2, 2008) 

The low energy structure of the energy momentum tensor in several chiral quark models with 
SU(3) flavor symmetry is analyzed in the presence of external gravitational fields. To this end 
a derivative expansion in a curved four dimensional space time is considered. This allows us to 
compute the chiral coefficients in, £12, and L13, appearing when the model is coupled to gravity 
and related to space-time curvature properties. 
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I. INTRODUCTION 

The energy momentum tensor plays a crucial role in 
quantum field theory since it arises as a Noether current 
of the Poincare group. It is conserved in all relativistic 
local theories even when there are no other conserved 
charges. Accordingto the standard Noether construction 
(see e.g. Refs. 0,12] for recent reviews on the classical 
theory), the energy momentum tensor can be computed 
as a functional derivative of the action with respect to 
an external space-time-dependent metric, g^ u {x), around 
the fiat space-time metric 77^ „ (we take the signature (-1 — 



)■ 



where 



s = 



ss 



-9 C{x) ■ 



(1) 



(2) 



The chiral Lagrangian C contains only metric contribu- 
tions, and it is obtained by promoting the flat space-time 
metric ij^^to the curved one g^ v . It takes the form given 
in Refs. H, 3- In QCD, the energy momentum tensor 
operator probes the interaction of quarks and gluons to 
gravitons. At the quantum level the high energy behavior 
of d^ v is improved if suitable additional transverse cor- 
rections are implemented |5(, and within such a set up 
there is a trace anomaly relating 8£ to the divergence of 
the dilatation current [y, 0, HJ , signaling the anomalous 
breaking of scale invariance. The non- vanishing vacuum 
expectation value, (O|0^|O) is related to the existence of 
a non- vanishing gluon condensate generating a bag con- 
stant and scale Ward identities QUE!- Deep inelastic 
scattering provides also some information on the momen- 
tum fraction carried by quarks and gluons in a hadron at 
a given scale fl2l | . Direct experimental determination on 



the basis of one graviton exchange is out of question due 
to the smallness of the gravitational constant as com- 
pared to weak and strong processes. The gravitational 
pion form factor can be used to determine the width of 
a light Higgs boson into two pions 01 > an d also the de- 
cay of a scalar glueball into two pions 0|. There have 
also been few attempts in the past to define on the 
lattice [l^ , but so far there are no results of practical in- 
terest concerning matrix elements between hadron states 
carrying different momentum. 

At low energies, the spontaneous breaking of chiral 
symmetry dominates and in the meson sector any opera- 
tor can be described in terms of a non-linearly transform- 
ing Pseudoscalar Goldstone boson field U with an infinite 
number of low energy constants (LEC's) |1 H 111 HllTst 
(for review see, e.g., Refs. 0,H3|). In a chiral expansion, 
the most general structure of 6„„ up to and including 
fourth order corrections reads jig 



9 - fl (0) 4- 6> (2) 

'flV — V I "pis 
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with 
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0(2) 
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D^D U U) - r,^C 



(2) 



(3) 

(4) 
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-r,^C^ + 2L 4 (D^D V U){ X ^U + U^ X ) 
L^DprfDJJ + DjflDJTitfU + U^x) 
2iii (r,^d 2 - dpd*) (D a U^D a U) 
2L 13 ( Vtw d 2 - d„d v ) ( X j U + U^x) 

x{D a U ] D f3 U) , 

(6) 



where (A) = trA means trace in flavor space and the 
expansion for the Lagrangian reads 
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£(0) + £(2) + £(4) + 



(7) 



The chiral expansion corresponds to an expansion in pow- 
ers of external momenta. The upper-script means the 
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order of the contribution. The zeroth order vacuum con- 
tribution reads 

£ (0) = B , (8) 

where B is the vacuum energy density constant. The 
second order contributions reads 

£ (2) = ^(DuUtLPU + tfU+rfx)), (9) 

and the fourth order contribution is 

= L^D^D^U) 2 +L 2 (D^D V U) 2 

+ L 3 {(D^D»U) 2 ) 

+ Li{D M UiD»V)(xiU + Uix) 
+ Ls(D^D»U(x ] U + U^x)) 
+ L^U + U^x) 2 

+ L 7 ( X ] U -U^xf +L S ((X^U) 2 + {U^x) 2 ) 
- iL 9 (Fj; v D> 1 UD v U' ! + F*D^U ] D V U) 
+ L W (F^UF^ R U^) 

+ ff!((^) 2 + (i^) 2 ) + ff 2 ( X t x ) . (10) 

Here, we have introduced the standard chiral covariant 
derivatives and gauge field strength tensors, 

D„U = D^U - U D* = d^U - i A^U + iU A*, (11) 
F; v = i[D^, Dl] = d^Al - d v Al - i[Al, Aft, 

with r = L,R. / is the pion weak decay constant in 
the chiral limit and U — e ,v ^ n /^ a unitary matrix and 
II a non-linear transforming field under the chiral group 
representing the octet of pseudoscalar mesons given by 









K+ 


7T 






> K ° 


K- 




K° 





(12) 



Finally, x = 2-Bo diag(m u , rrid, nris) is the current quark 
mass matrix with Bq — \(qq)\/f 2 where (qq) is the quark 
condensate for one flavor. 

Note that the coefficients L\ — Liq appear in given 
by Eq. IjlUfl . As one sees, the terms containing L\\ — L13 
cannot be obtained by computing the energy momentum 
tensor from the chiral effective Lagrangian in flat-space 
time and from this viewpoint are genuine new contribu- 
tions to 9^". The gravitational low energy coefficients 
arise at the hadronic level due to quantum effects and 
corresponds to the inclusion of curvature terms in the La- 
grangian which generate purely transverse terms in the 
energy momentum tensor. Actually, from a calculational 
viewpoint it is more advantageous to couple the system 
to gravity and compute the chiral Lagrangian in curved 
space-time. 

Although the LEC's should be deduced from QCD it- 
self solely in terms of Aqcd and the quark masses, their 



numerical values are obtained from a fit to low energy 
data, and systematic calculations of the corresponding 
non-gravitational low-energy constants (LEC's) Li_io 
have been carried out in the recent past up to two loop ac- 
curacy I2lll2^l2p . l2^ |. or by using the Roy equations [2{| 
(see also |2a|23|). For strong and electroweak processes 
involving pseudoscalar mesons the bulk of the LEC's is 
saturated in terms of resonance exchanges |2*sll2*9ll3cLl3l| . 
which can be justified in the large- N c limit in a certain 
low-energy approximation [82l | by imposing relevant QCD 
short-distance constraints. In the case of gravitational 
processes similar ideas apply ^8j, although less model 
independent information is know n. ( A one loop calcula- 
tion has been carried out in Ref. 123). 

In a more model dependent set up the LEC's can also 
be understood from a more microscopic view p oin t by us- 
ing chiral quark models |3^ . l35t l36t I37L l3St 1391 |40L l4lL 14^ . 
13 EE EE E3 for the non-gravitational coefficients 
L\ — Liq. The calculation of Ln,Li2 and L13, encoding 
the coupling to gravitational sources, has only been done 
for a chiral quark model in Ref. |48j and more recently 
in the spectral quark model |49j. These values could be 
used to predict the experimentally unmeasured gravita- 
tional form factors of pseudoscalar mesons (see Ref. |3^| 
for a calculation within ChPT). 

In this paper we provide the values of the gravitational 
LEC's £11,12,13, as well as i?o,3-5, in several chiral quark 
models. Namely, we use in Sect. [H] the Nambu-Jona- 
Lasinio (NJL) model [Hoi l5ll | and the Georgi-Manohar 
(GM) model |52j in the presence of gravity. Although 
these models look very different in appearance they do 
incorporate chiral symmetry breaking and quarks have 
a constituent mass M ~ 300MeV. The main difference 
between them is that while in the NJL model this mass is 
dynamically generated through interactions, the GM al- 
ready assumes that this breaking has taken place. More- 
over, while in the GM model there are only pseudoscalar 
Goldstone bosons as dynamical degrees of freedom in the 
NJL there are additional qq scalar fields. Finally, in the 
NJL model quarks appear to have an axial coupling con- 
stant qa for the quarks equal to one, while in the GM 
model one starts already with qa 5^ 1. A more de- 
tailed analysis ^1 shows that if in NJL Gy ^ then 
gA 7^ 1 (see also below). More specifically, we consider 
the generalized SU(3) NJL with scalar, pseudoscalar, vec- 
tor and axial couplings extending the original calculation 
of Ref. [42j (see also Ref. [45j where also an estimate on 
the gluonic effects of the LEC's was done.). To this end 
we apply, after bosonization and Pauli-Villars regular- 
ization, a coordinate and frame invariant derivative ex- 
pansion by means of the heat kernel technique in curved 
space-time (Sect. After that, the scalars, vectors 

and axial vectors are eliminated at the mean field level 
fSect. ITVjl and the full results are presented in Sect. El 

For notation details, we rely heavily on previous 
work 0j 1491 E3|. Nevertheless, the inclusion of grav- 
ity is a bit messy, so we present the minimal amount of 
formulas to make the paper more comprehensive and also 
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to point out some differences with previous calculations 
non-gravitational LEC's by one of us (ERA) [4^1 and 
other authors as well as our previous work with Bro- 
niowski on the gravitational LEC's in the spectral quark 
model Hij. 



II. CHIRAL QUARK MODELS IN CURVED 
SPACE TIME 

The coupling of gravity 1 to chiral quark models has 
been described at length in Ref. ^Ij where more de- 
tails on the notation and this specific model can be 
found in the context of the Spectral Quark Model 
(SQM) |5^, |5|| In this paper we discuss several 

constituent chiral quark models which have the com- 
mon feature of incorporating dynamical chiral symmetry 
breaking at the one quark loop level: The Nambu- Jona- 
Lasinio (NJL) model [5(1 l5lj | and the Georgi-Manohar 
(GM) model [5i| • In these models quarks have a con- 
stituent mass M ~ 300MeV. The main difference be- 
tween them has to do with the presence or not of dy- 
namical qq scalar fields respectively. In addition, while 
the NJL model dynamically generates spontaneous chi- 
ral symmetry breaking, the GM model already starts in a 
chirally broken phase. To avoid unnecessary duplication 
we will discuss explicitly the NJL model with scalar and 
vector couplings from which the relevant implications for 
the GM model can be more easily deduced. 



A. Nambu— Jona-Lasinio model 

The Nambu- Jona-Lasinio (NJL) action in curved 
Minkowski space-time with metric tensor g^ u (x) reads 



S = J d xy/^gC N j L (13) 
where g = Det(<7 At!y ) and the Lagrangian is given by 



^ £ ((qXal.q) 2 + {qXaWqf) (14) 



where q = (u,d, s, . . .) represents a quark spinor with 
N c colors and Nf flavors. The A's represent the 



1 We consider only Einstein gravity, i.e., we use the Riemann con- 
nection, uniquely defined by being torsionless and metric pre- 
serving (metricity). Extension to torsion gravity l- r )4l is possible 
but will not be studied here. 



Gell-Mann matrices of the U(Nf) group and Mq = 
diag(m u , ma, m s , . . .) stands for the current quark mass 
matrix. In the limiting case of vanishing current quark 
masses the classical NJL-action is invariant under the 
global U(Nf)n (g) U(Nf)L group of transformations. The 
derivative — iu^ is frame (local Lorentz) and general- 
coordinate covariant and it includes the spin connection, 



■r 



(15) 



where the frame covariant derivative is defined as usual 



7«iP = d^ v {x) - r v Ax)l\[x) 



(16) 



and are the Christoffel symbols which for a Rieman- 
nian space read, 



1 



r V = i d ^9^ + d^g\ v - d v g^} 



(17) 



The space-time dependent Dirac matrices, 7 M (x), fulfill 
the standard anticommutation rules {7^(2;), 7^(2;)} = 
2g flu (x). The 75 matrix is, on the contrary, x- 
independent (we use the conventions of Ref. 58] ). The 
Dirac slash is constructed as i</Dq{x) — ij fJ, (x)d t _ l q(x) 2 . 

The vacuum to vacuum transition amplitude in the 
presence of external scalar, pseudoscalar, vector, axial 
and metric bosonic fields [s,p, v, a, g) Lagrangian can be 
written as a path integral in the form 



Z[s,p,v,a,g] = (0|Texp|i J d 4 Xy/^g£ cxt (x)} |0) 



DqDq exp|i 



NJL 



where the external field Lagrangian is given by 
£ cxt (a;) = q(j + ^75 - (s + ijap)) q ■ 



- £ext] j , 
(18) 

(19) 



Following the standard procedure |59j| it is convenient 
to introduce auxiliary bosonic fields (S, P, V, A) so that 
after formally integrating the quarks out one gets the 
equivalent generating functional 



Z[s,p,v,a,g] 



DSDPDVDAe 



ir[S,P,V,A;s,p,v,a, : 



(20) 



2 Note that in the curved case d^ip(x) = g>* u (x)d v <p{x), and the 
self-adjoint D'Alambertian on a spin-0 field is given by d 2 ip(x) = 
l/^f—gd^^f—gg^dv^pix)) so that 

d A x^J~^g~ip\ (x)d 2 ip2(x) = J d 4, xyj~^gip2(x)d 2 ^p\{x) . 

Likewise, for a Dirac spinor the form J d A X\/— gi/>i(i^ + i/>)t/>2 is 
Hermitian. 
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where the effective action 

r[S,P,V,A;s,p,v,a,g} = -iV c Tr log (zD) 



(21) 



have been introduced 3 . Here, the Dirac operator is given 
by 

iD = i@ + y!,-M + (f + 475 -S-i^P) 
where we introduce the additive combinations 



S = s + S, P = p + P, V = v + V, A = a + A 



(22) 



/(A 2 ) + A 2 /'(A 2 ). In the vacuum (s=p = v = a = 0) 
one has spontaneous chiral symmetry breaking, since by 
minimizing r cvcn one gets S = M with M the constituent 
quark mass and P = V = A = 0, and the gap equation 



1 

Gs 



x - ; d 4 k 
= -iAN c 2_^c 



1 



(2tt) 4 k 2 - M 2 - Af 



(26) 



where the Pauli-Villars regularization has been used. 
Thus, according to Goldstone's theorem we will write 
the decomposition 

M = S + iP = MU + m (27) 
with U = e lv/2n / / and II is given by Eq. Q . 



and the mass term 



r — 



4G. 



-(S 2 + P 2 ) 



4G, 



.(V^ + A fl A^)(23) 



has been defined. Here, (5, P, V, A) are dynamical scalar, 
pseudoscalar, vector and axial meson fields respectively 
whereas (s,p,v,a,g) represent external sources. With 
the exception of the metric tensor g^(x), all fields are 
expanded in terms of the A flavor matrices. Notice also 
that for the path integral in the bosonic fields to be well 
defined in Minkowski space we must use the prescription 
Gg 1 -» Gg 1 - ie and G v x -» G v x - ie. 

The additive contribution to the effective action can 
be separated into a 75-even and 75-odd parts correspond- 
ing to normal and abnormal parity processes respectively. 
Introducing the operator 4 



B 5 [S,P,V 7 A}^ J5 B[S,~P,V,-A} 75 



(24) 



the 75-even contribution is quadratically divergent and 
can be regularized in a chiral gauge invariant manner by 
means of the Pauli-Villars scheme 

Tcven = -^Tr^ Cl log(DD 5 +A 2 +ie) (25) 

~4(b (52 + ^ + 4G7 ( ^ + 

where the Pauli-Villars regulators fulfill cq = 1, Ao = 
and the conditions J+ c$ — 0, J+ qA 2 = which render 
finite the logarithmic and quadratic divergences respec- 
tively. In practice, we take two cut-offs in the coincidence 
limit Ai — > A 2 = A and hence X« C »/( A ?) = /(°) _ 



3 For a bilocal (Dirac and flavor matrix valued) operator A(x, x 1 ) 
one has 

j4„ 



Tr A = / drx-y t /-gti{A(x, x)) , 

with tr denoting the Dirac trace and ( } the flavor trace. 
4 D5 corresponds to rotate D to Euclidean space, take the hermi- 
tian conjugate and rotate back to Minkowski space. 



B. Georgi-Manohar model 

In the presence of gravity the Georgi-Manohar model 
Lagrangian j^l reads 



£gm = q[ i@ + 0-MU 5 -M Q 



+ ^(l-gA)U 5 i0U 5 \q 



(28) 



where U 5 = e i(v / 2n//) 75 ) and n is given by Eq U^J), M 
the constituent quark mass, and gA the axial coupling 
of th e q uarks, which will be assumed as suggested in 
Ref. [52] to be different from one. The generating func- 
tional now reads, 



Z[s,p,v,a,g] = J DUe lT[u '' s ^ v ' a ' 9] 



(29) 



with the effective action given by 

T[U;s,p,v,a,g] = -LY c Trlog (»D) . (30) 

Direct comparison with the one quark loop effective ac- 
tion of the NJL, Eq. 121|) shows that it corresponds to a 
similar model without mass terms C m and with a Dirac 
operator as in Eq. (1221) with a specific choice of the spin-1 
dynamical field operators, 

V, = ~(1 - 9a) [U^d.U - d^U] , (31) 
A^ = ^(l-g A )[U^d f2 U + d^U] . (32) 

In Eq. (|3l)[l we understand that the same Pauli-Villars 
regularization as in the NJL model has been imple- 
mented. 
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III. HEAT KERNEL EXPANSION 

In the chiral expansion of the action, Eq. i|21[l. the 
counting is the standard one, the pseudoscalar field U 
and the curved space-time metric are zeroth order, 
the vector and axial fields and a M are first order, and 
any derivative is taken to be first order. The external 
scalar and pseudoscalar fields s and p and the current 
mass matrix mo are taken to be second order. At the 
one quark loop level this chiral expansion corresponds to 
a derivative expansion which should be gauge, frame and 
coordinate invariant. 

In order to carry out the low energy expansion of the 
action via a heat kernel method it proves useful to use 
the proper-time integral representation of the logarithm 



E 

i 



r°° ri-r 

~] dTr log (D 5 D + A?) = -Tr / -e"- 13 ^ 

; JO T 



where <P(t) is the proper time representation of the Pauli- 
Villars regularization, 



E 



cue 



with the conditions 

5> = o 



E c ^ 2 = o 



(34) 



(35) 



which fulfills the conditions 0(0) = and </>'(0) = 0, thus 
killing the quadratic and logarithmic divergencies respec- 
tively. The operator inside the logarithm is of Klein- 
Gordon type in curved space-time with some spinorial 
structure, 



D,D 



f 



[VvW^gsTVu)] +v, 



with 



and 



V = VrPr + V L P L 



v^F* + 7 R- h^^M + M^M, 



(36) 



(37) 



(38) 



As we saw in the previous section, the one quark loop con- 
tribution to the effective action depends additively on the 
external and dynamical fields. Although this complicates 
matters one can borrow from previous results taking into 
account some minor modifications regarding this fact and 
the regularization 5 . 



5 This additivity works only for the normal parity contribution of 
the action, but must be modified in the abnormal pari ty secto r 
in order to reproduce the QCD anomaly. See e.g. Ref. [53, 61] 



The form of the operator D5D in Eq. (|36() is suitable to 
make a heat kernel expansion in curved space-time as the 
one of Ref. pSj. F° r a review see e.g. and references 
therein. In our particular case, before undertaking the 
heat kernel expansion we separate a constituent quark 
mass squared, M 2 , contribution from the operator D5D 
which we treat exactly, 



4tD,D I 



rM 2 



-ir(D 5 D-M 2 )| 



(39) 



-irM 



(47rir) 



n=0 



The r integrals appearing in the effective action are of 
the form 



T21 ■= M 



dr 

— 0(r)(ir)'e 



l -irM 2 



(40) 



For completeness, we list here the particular values, 
M 4 Z_4 = E C *( A " + M ') 2 lQ S( A2 + M ') > ( 41 ) 

i 

AI 2 l- 2 = Yl + A/2 ) lo §( A2 + A/2 ) > ( 42 ) 

i 

X = -^c,:log(A 2 + Af 2 ), (43) 

i 

M 2 



2 2n = r(n)2j. 



A 2 + M 2 



, Re(ra) > (44) 



After evaluation of the Dirac traces, the second order 
Lagrangian is given by 

(4tt) 2 

M 2 
~6~ 

whereas the fourth order becomes 

(4) = _Nc_r 

9 (4tt) 2 I 



4 2) = ^{ M2 M%U^U) (45) 

M 2 1 
+ 2M 3 I_ 2 (m t [/ + U^m) + —1- 2 (R) \, 

6 J 



(4tt) 2 

- \M{F^) 2 + {F^) 2 
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1 



1 



l 720 



144 90 



1_ „=« 



^Z 4 ((V M £/V,^) 2 ) - ^((V^V'V) 2 ) 
+ 2M 2 X_ 2 (7n t m) - M 2 l ((fn''U + U^m) 2 ) 

- MX 2 (v M ^ t v M L/(7fr t c/ + c/ t 7fr)) 

+ MIo (V M U^m + V p m f U) 

- ^-X R{tfm + M*U) - ^R^U^U)}. 

(46) 
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Here R^^, Rfw and R are the Riemann curvature tensor, the following relations, 



the Ricci tensor, and the curvature scalar, respectively , 



pA pi -p\ o -pA 



R/iu — R 



t^A pa _ -pA "pet 
llol va vol fia") 



R = g 



[IIS ■ 



(47) 



The Riemann connection is given by the Christoffel sym- 
bols, defined in Eq. (|T7|l. 

We have introduced the gauge and frame covariant 
derivatives and the field strength tensor containing the 
external and internal (bosonized) degrees of freedom, 



V„17 



iVfU + iUV? 



(48) 



with r = L,R and the spin additive combination 

m = (S + iP- MU) + m . (49) 

The form of Eq. i|46|l is not yet ready for comparison with 
the result of Ref. [4J, [l8( . To do that we have to eliminate 
all degrees of freedom other than the on-shell pions. We 
proceed in three steps: Firstly we integrate out the vector 
and axial degrees of freedom, afterwards we eliminate the 
scalars and finally we exploit the classical equations of 
motion for the pseudoscalars. 

IV. MEAN FIELD EQUATIONS OF MOTION 
A. Elimination of Vector and Axial Vectors 



To eliminate the and A^ fields at the mean field 
level we minimize the corresponding Lagrangian with re- 
spect to those fields. To the order in the chiral expansion 
we are working it is enough to deal with the terms con- 
taining vector mesons with two Lorentz indices, i.e. the 
mass term and the second order term stemming from the 
quark determinant, 

4 2 V = j^M'loiV^V^ + ^-iV^ + A^) 



(4tt) 2 



'V 



F R 

(XI/ 




-9a)FK+\{1- 






4 l 


- a 2 *) (V ',,U^V„U 

•j A. 1 V M v 


-V„U^V„U) (53) 


F L 

(11/ 


-!<>- 


-g A )UF*tf + \ 






i , 

~ 4^- 


-<&)(v /1 tfv„i7 t 


-VuUVJJ*) ,(54) 


V M /7 


= 5aV p 


u, 


(55) 



VV = g A V 2 U + ig A (l-g A )UV^WU. (56) 



B. Elimination of Scalars 

The elimination of scalars proceeds along similar lines 
as in the vector and axial case. If we do a chiral rotation 



s + iP = VuzVu 



(57) 



with £t = S, and using that E = M + $ with $ the 
fluctuation around the vacuum value we have 



m = Vu^Vu + — =-x ■ 
2£>o 



The mass term becomes 

C - 1 



4G S 



(M 2 + 2M$ + $ 2 > 



(58) 



(59) 



(50) 



The result looks like the GM model Eqs. (|31I32|) as noted 
in Ref. H2, 



v R 



V, 



v R + l -(l-g A )U^^U, 



vfl + -(l-g A )UV^ 



(51) 
(52) 



which looks like a chiral gauge transformation of the ad- 
ditive fields V and A. Then it is straightforward to obtain 



Note the opposite si gn o f our definition f or t he Riemann tensor 
as compared to Ref. lid . We follow Ref. |Q4| . 



Using the gap equation, Eq. (|26[1 . linear terms in <f> not 
containing external fields vanish. As a consequence, the 
part of the Lagrangian comprising the scalar field is given 
by 

U(x) = --^(^M 2 X ^ 2 + ^MX R^ 

+ MioVu^Vw (uW 2 u f + v 2 uuA 

+ 2M 2 (2Z - X_ 2 )Vc7$Vc7t( I 7m t + U^m) 
+ 2MZ 2 VZ7$v / CnV AI [/V' 1 [/ t ^ (60) 

from which the well known mass formula Ms — 2M fol- 
lows. Minimizing with respect to the field $ the classical 
equation of motion follow, 

Substituting this equation into the Lagrangian we 
obtain the contribution to the effective Lagrangian stem- 
ming from the integration of scalars. 
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C. Mean field equations for pseudoscalars 



and 



The relevant equations of motion for the non-linear U 
field are obtained by minimizing C^ 2 \ One obtains a set 
of relations which are valid in the presence of curvature, 

(v 2 c/tv 2 t/) = ((v^t/Vc/) 2 ) - i ( ( x t [/ _ c/ t x ) 2 ) 



(62) 



and 



( x W 2 u + v 2 u^ x ) = ^x)-\{{x^u + u\f) 

- ((x t C/ + C/ t x)V^tv M C/) 
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(63) 



C {i ' 9) = LiiV^VU) 2 + L 2 (V A1 C/ t V„[/) 2 

+ L 3 ((V,U^»U) 2 ) 

+ uiw^w^ix^u + u^x) 

+ L 6 (x f U + U^ X ) 2 

+ LrtfU-rtxf+LsitfUf + irtx) 2 ) 



+ L W (F^UF^ H W) 



(69) 



In the case of the C/(3) group one has Det U = e iri0 ^ 
not necessarily equal to unity and the last two terms 
involving (x^U ± U^ X ) 2 in Eqs. JHU an( i <E3 should be 
dropped. There is another integral identity which proves 
very useful 

J tfW^V^V^V^Vt/) = J d i x^{S/ 2 U i V 2 U) 
d 4 x^(F^UF^ R U^) 



(64) 
(65) 

2 



The curvature contributions to the chiral Lagrangian can 
be written in the form proposed in Ref. [Tsl ] and are given 
by 



= -HnR 



(70) 



and 



£ (4,H) 



U + U^x) + H3R 2 + HiR^R^ 
r.-apR^. (71) 



-^13 

H 5 R 



Finally, we also have the SU(3) identity 
((V^^V,?/) 2 ) = -2((V M f/ t V^C/) 2 ) 



+ (V M C/ t V„C/) 2 + -(V M C/ t V^C/> 



The curvature terms reflect the composite nature of the 
pseudoscalar fields, since in the considered model they 
correspond to the coupling of the gravitational external 
field at the quark level. 

Note that the pseudoscalar matrix U is a coordinate 
and frame scalar. So, once and only once the identities 
l|55 )l -l(B3 )l have been used one can substitute coordinate- 
frame covariant derivative by the covariant derivative, 
i.e., V M £7 = Dfj,U. In this way Eqs. © and ljTU|) of this 
manuscript, and Eq. (10) of Ref. [l8j are deduced. 



GM model 



V. RESULTS FOR THE 
GASSER-LEUTWYLER-DONOGHUE 
COEFFICIENTS 

Following the steps of Ref. one S e ^ s the chiral ef- 
fective Lagrangian in the presence of gravity in the form 
of Gasser, Leutwyler and Donoghue form p.ll^| 

The zeroth order vacuum contribution reads 



£(0) = B = _2NfNc Mil _ 
(Airy 



(67) 



where B is the vacuum energy density constant. The 
metric contributions read 



£(2, 5 ) 



r- 



(V^V^U + ^U + U^x)), (68) 



The Gasser-Leutwyler-Donoghue coefficients for the 
GM model do not contain any contributions from the 
scalars, i.e. the spin field. So, the only contribution 
stems from the quark loop. For this model, the pion 
weak decay constant is 



/ 2 = ^1M%. 
The normalization factor for the field x is 

M l- 2 



B 



9a j o 



With 



M f 2 , 

B I (qq) I 1-2 



(72) 



(73) 



(74) 
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the result we find for the GLD coefficients reads 



u = 



48(4tt) 2 
+ 4g 2 A (l - g 2 A )X 2 + 2g 4 A J i \ , 
Li = 2L X , 

AT 

3(l-g 2 A ) 2 l + 8g A l 4 



L 3 

Li 
L 5 
L 6 
L 7 

Ls 

Lg 

Lio 

Ln 

L\2 
L13 

Ho 
Hi 
H 2 
H 3 
Hi 
Hn 



N r . 



24(4tt) 2 



i(3 - 4ffi)X 2 



0, 



Nr 



2(4tt) 2 
0, 



P9a [La - L 2 \ 



Nr 



24(4 7 r) 2 iV / 
N 



g A [6pl - g A L 2 ] 



24(4tt) 2 
N 



[6p(p - g A )J + g A L 2 ] 



6(4tt) 2 
N c 



6(4vr; 
N c 



f2(47r) : 



[(1 - g 2 A )I + 2g 2 A I 2 ] . 

-[(l-g A )l + g\l 2 ] 



2 T 



N 



6(4tt) 2 
N 



9A J 2 , 



f2(47r) : 

N r N 



;pLo 



P P 

ASM 2 g\ ' 
Nff 2 



M z I- 2 = J -—, 

24 p ' 

-(l+g 2 A )l +g 2 A L 2 ] 



6(4tt) 2 

N c 
12(4tt) 2 

N c 
12(4tt) 2 
N c N f 

f44(47r) : 
NcN f N f f 2 



(75) 
(76) 

(77) 
(78) 
(79) 
(80) 
(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 



[6p 2 X_ 2 - 6p(p + g A )l + g\L 2 \ (90) 



Nf P 
576M 2 g\ ' 



90(4tt) 



360A/ 2 g\ 



lN c Nf = 7N f f 2 
720(4tt) 2 2880M 2 <? 2 1 



(91) 
(92) 
(93) 



We take M = 300 MeV and g A = 0.75. For given values 
of M and g A , the cutoff A is adjusted to reproduce f n = 
93.2 MeV. This yields 

A = f 470 MeV , B = 4913 MeV , 

I- 2 = 20.8 , J = 2.26 , 1 2 = 0.922 , 

J 4 = 0.995 . (94) 

The constituent chiral quark model (QC) corresponds to 
take g A = 1 in the previous coefficients. Using the same 



value for M, we have for this model 

A = 828 MeV , B = 1299 MeV , 

Z_ 2 = 5.50, T = 1.27, 22 = 0.781, 

Ji = 0.963 . (95) 

The numerical values for the GLD coefficients are dis- 
played in Tabic [I] 



B. NJL model 

The GLD coefficients for this model will have two dif- 
ferent contributions: one coming from the quark loop 
and subsequent integration of spin f fields, and another 
coming from the integration of spin fields. For the 
quark loop contribution we have the same expressions as 
Eq. H75|) - (|93|l . The pion weak decay constant is 



1 N r o 

f = ^9aM 2 1 q . 



(96) 



Note the power g A as opposed to the power g 2 A in the 
GM model, Eq. {721 . T ^ c difference has to do with the 
absence of a mass term like C m in the GM model. We 
use the notation 



Bo 
With 



M 



2G s p 



= -^, g A = l-2G v f 2 . (97) 
9a L 



M 
B~o 



9A- 



the spin contribution becomes 
N, 



Lo = 



4(4tt) 2 To 
N 



9a [Lo-L 2 ] 2 , 



2 9a(9a 



^8 



^ll 



-^13 

Hi 



H> = 



4(4tt 

N c 
16(4tt) 

N c 
f2(47r) 2 

N c , 
2^f {9A 

N c Nf 
f44(47r) 



2p)[lo-L 2 ] 
2 (.9 a - 2p) 2 I , 



9a [Lo - L 2 \ 



Nf P 
576M 2 g A 



(98) 

(99) 

(100) 

(101) 

(102) 

(103) 
(104) 
(105) 



The remaining Lf , Hf vanish identically. The sum of 
both contributions will give the GLD coefficients for the 
NJL model, which read 



Li = 



N r 



48(4tt) 5 



(1 



+ 4 5 2 (f- 5 2 )Z 2 + 2< ? lZ 4 



(106) 
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U 

U 
L 5 
L 6 
L- 



1L X 



N c 



24(4tt) 2 
2(3- 



3(1 - 2g\ - g\)I + : 



3a) -3^— )1 2 



0, 



Nr. 



4(4tt) 2 
0, 



g\[l -l 2 ] , 



24(4tt) 2 V / 



9A [6p2" - g A l 2 ] 



L H = 



La = 



Lio — 



N n 



48(4tt) 2 

N c 
6(4tt) 2 
N c 
6(4tt) 2 
N 



g\ [3X - 2Z 2 
[(1 - g 2 A )X + 



[(l-g 2 A )l +g 2 A l 2 ] 



12(4tt) 2 



6(4tt) 2 
N 



9a t o z 
9a j 2 ■ 



-^13 = 

Ho = 

H x = 

H 2 = 

H 3 
Hi 



24(4?r) 2 
N r N 



gA?o 



9aP 
ASM 2 



P 
96M 2 



6(4tt) 2 

N c 
12(4tt) 2 

N c 
24(4tt) 2 



f 9 N f p 



24 p ' 
[-(l + g 2 A )l +g A l 2 ] , 

[l2p 2 X_ 2 + 3g A (g A - 8p)l Q 



0, 

N c N f 
90(4tt) 2 

7N c N f 
720(4tt) 



Xn = 



360M 2 g A ' 
7N f P 
2880A/ 2 g A 



(107) 

1 

(108) 
(109) 
(110) 
(111) 
(112) 

(113) 

(114) 

(115) 

(116) 

(117) 

(118) 

(119) 

(120) 

(121) 
(122) 

(123) 
(124) 



The coefficients L\, L 2 , L4, Lq, L7, Lg, L X q, L\ 2 , Hq, 
Hi, H4 and H$ in the GM model coincide with those 
of the NJL model. However, we prefer to display them 
explicitly for easier reference and because the expressions 
for p (Eqs. O and @) do not coincide. 

Note that this model reproduces the relation L3 = 
— 6L1, provided terms 0(N c g A ) are neglected. We note 
several differences with previous works. The values L\, 
L 2 , L3, L4, Lg, Lq, Lg, Lio, Hi and H 2 coincide with 
Ref. |45| . Lg, differs in two powers of g A in the term 
proportional to T 2 . (We reproduce their results for ev- 
ery separate pieces: quark loop contribution and spin 
contribution.) 

The value of L7 is non-cero, if one imposes the condi- 
tion Det(U) = 1 within SU(Nf) flavor symmetry. Both 
in Ref. |42j and |45| this term is not obtained, in spite 



of the fact that in those papers it is explicitly mentioned 
that one works in SU(iV^) flavor symmetry. The situation 
can be mended if one considers U(iV/) flavor symmetry 
instead where Lj = 0. 

The coefficients Li_io where given in Ref. [42]. The 
present values of L4, L5, Lq, Lg,, Lg and Lio coincide with 
those |42| where an extra erroneous term in L\ appears. 
L3 differs from Ref. 42] in all factors except for the one 
in X4. We correct here these errors (see also Ref. |45|). 
Hi and H 2 did not appear in that reference. 

The coefficients Lu, L i2 and L13, as well as H , 3 ^ 5 , 
are new and are the main result of this work. They have 
been also evaluated some time ago ^8| in a quark model 
without scalars and vectors and more recently by us 0] 
in the spectral quark model. 

The numerical values for these coefficients Eq. <|106f) - 
(|124f> are displayed in Table U in two different cases: one 
for the generalized SU(3) NJL model, and another one 
in which integration of spin 1 fields is not considered, 
i.e. g A = 1. For the first case we take g A — 0.606 as a 
reasonable value. Using M = 300 MeV, this yields 



A = 1344 MeV , B = 4015 MeV , 
Z_ 2 = 17.0, T = 2.10, 1 2 = 0.907, 
X 4 = 0.993 . 



(125) 



On the other hand, for g A — 1, the numerical values of 
A, Bo, p and X 2n are identical to those for the QC model, 
Eq.©. The LEC's in the NJL model with g A = 1 and 
in the QC model differ due to the scalar contributions 
Lf 58111 3 and ff 2 s 3 which are not present in the QC 
case. 



VI. SUMMARY AND CONCLUSIONS 

In the present work we have calculated the low energy 
constants of the energy momentum tensor. Some of these 
LEC's are directly determined by the standard Gasser- 
Leutwyler coefficients, whereas other, Ln_i3 and iLj,3-5j 
are new as are driven by operators not present in the (flat 
space) chiral Lagrangian. Technically, the best way to 
proceed is to consider QCD in a curved space-time be- 
cause these allows us to work with a single scalar, the low 
energy Lagrangian, rather than its variation, the energy 
momentum tensor, making easier both the computation 
and the imposition of symmetry restrictions. In addition 
the coupling to gravity is also obtained as a byproduct. 
The curved space chiral Lagrangian, contains two kinds 
of contributions. On the one hand, those obtained by 
"minimal coupling" from the flat space case, £^ 9 ' , on the 
other, non-minimal pieces containing the Riemann cur- 
vature tensor C^ R \ In the spirit of not introducing new 
fields other than the metric tensor, we have considered 
only Einstein's gravity. Further new pieces would appear 
by using connections with torsion or violating metricity. 
As happens with gauge couplings (e.g. the magnetic mo- 
ments) the non-minimal gravitational pieces are not fixed 
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TABLE I: The dimensionless low energy constants and -ffo compared with some reference values and other models. The values 



quoted for Li_io,-Hi_5 are 


to be multiplied by 10 3 . 


The value quoted for Hq 


is to be multiplied by 10 3 MeV 


2 




ChP 1 


NJL 


NJL 

f 1 \ 

(9 A = 1) 


QC 


GM 


(MDM) 


Large N c c 


Dual 6 
Large N c 


Li 


0.53 ± 0.25 


0.77 


0.76 


0.76 


0.78 


0.79 


0.9 


0.79 


L> 


0.71 ± 0.27 


1.54 


1.52 


1.52 


1.56 


1.58 


1.8 


1.58 


L z 


—2.72 ± 1.12 


—4.02 


—2.73 


—3.62 


—4.25 


—3.17 


—4.3 


— ^ 1 7 


La 


























Lr, 


0.91 ± 0.15 


1.26 


2.32 


1.08 


0.44 


2.0 ± 0.1 


2.1 


3 1 7 


L 6 























n 

u 


L- 


—0.32 ± 0.15 


—0.06 


—0.26 


—0.26 


—0.03 - 


0.07 ± 0.01 


—0.3 




L s 


0.62 ± 0.20 


0.65 


0.89 


0.46 


0.04 


0.08 ± 0.04 


0.8 


1.18 


Lg 


5.93 ± 0.43 


6.31 


4.95 


4.95 


6.41 


6.33 


7.1 


6.33 


Lio 


—4.40 ± 0.70 d 


—5.25 


—2.47 


—2.47 


—4.77 


—3.17 


—5.4 


—4.75 


in 


1.85 ± 0.90 c 


1.22 


2.01 


1.24 


0.82 


1.58 


1.6 e 




il2 


— 2.7 e 


— 1.06 


—2.47 


—2.47 


— 1.64 


—3.17 


— 2.7 C 




-^13 


1 . 1 zn u.ou 


1 01 


1 m 


47 


0.22 




l.l e 




H 




-14.6 


-4.67 


-4.67 


-17.7 


1.09 






ffl 




-4.01 


-2.78 


-2.78 


-4.76 








H 2 




1.46 


1.45 


0.59 


0.49 


-1.0 ±0.2 






H 3 










-0.50 


-0.89 








H 4 




1.33 


0.80 


0.80 


1.43 








H 5 




1.16 


0.70 


0.70 


1.25 









"The two-loop calculation of Ref. |2 3| 



^Ref. 


m 


c Ref. 


M 


d Ref. 


M 


e Ref. 


m 



(nor demanded) by general covariance of the chiral La- 
grangian and to obtain them requires to couple gravity 
directly to QCD quarks and gluons previously to their 
integration out to yield the low energy Lagrangian. 

The standard LEC's, which are expected to follow from 
QCD, have also been computed in QCD-like models con- 
taining chiral quarks, with some degree of success. Here 
we apply the same approach for the non-minimal pieces. 
We do so for various models, namely, the constituent 
quark model, the NJL model with and without vector 
mesons, and the Georgi-Manohar model. This extends 
our previous calculation carried out within the spectral 
quark model 0] where also a dual large N c model was 
introduced. The results quoted in Table [I] include stan- 
dard and non minimal coefficients in these models as well 
as ChPT calculations and the single resonance Large N c 
model. The results for the LEC's look roughly quite simi- 
lar. As a rule all models and fits give the same sign for all 
coefficients, with the exception of Hq and H2 in the SQM. 
For the standard Gasser-Leutwyler coefficients Li_io the 
best overall agreement with the two-loop ChPT calcula- 
tion of is achieved by the NJL model with vector 
mesons, for which the value of the reduced chi-square 
X 2 /N, (N = 10), is 2.2, although the QC and GM mod- 
els give results of near quality: 2.5 and 3.6 respectively. 



For the new coefficients there are no widely accepted val- 
ues. The closest agreement with the resonance satura- 
tion and large N c estimates of for _Ln_i3 is provided 
by the NJL model without vectors mesons, for which 
X 2 /3 = 0.29, but this is, of course, not conclusive. We 
should also mention the remarkable agreement between 
the prediction of the SQM for these three coefficients and 
that of the chiral bosonization model of . 



APPENDIX: A SPECTRAL QUARK MODEL 
WITH g A 1 ? 

One of the common features of both the NJL model 
with Gy and the GM model analyzed in this pa- 
per is the presence of an axial coupling constant for the 
quarks gA 7^ 1- In this appendix we analyze the possi- 
bility of extending the Spectral Quark Model to have 
a coupling similar to the GM model. Let us remind 
that the SQM [5^, |5^, is characterized by a gener- 
alized Lehmann representation for the quark propagator 
in terms of a spectral function p(uj) defined on a suitable 
contour on the complex mass plane, which according to 
the rules devised in those works should satisfy the spec- 
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tral conditions for the positive moments 



dusp(w)u n — 6 n Q , n > . 



(A.l) 



The vanishing of even moments precludes a positive def- 
inite spectral function. These conditions arise because 
one demands that high energy processes are free from 
logarithmic corrections as it corresponds to an effective 
model at low energies. This is the case for the electro- 
magnetic form factor of the pion or the partonic distri- 
butions in the pion. Gluonic radiative contributions can 
be implemented by renormalization group equations us- 
ing the quark model results as initial conditions for the 
QCD evolution (see e.g. Ref. |6t|h 

The physical hadronic observables turn out to depend 
both on the negative moments and the so-called log- 
moments, 



(A.2) 



The set of spectral conditions, Eq. (jA.lfl . remove the de- 
pendence on the scale /i in Eq. l|A.2fl and thus prevents 
the occurrence of dimensional transmutation for all mo- 
ments except p' . Assuming that the quark model is de- 
fined on a very definite renormalization scale, any low 
energy hadronic observable should not depend on any 
additional arbitrary scale, so we expect that p' does not 
determine hadronic properties, and in particular low en- 
ergy constants corresponding to process involving pseu- 
doscalar mesons, i.e. L\ — £1 3. Actually, in our low 
energy analysis of the SQM [4!| the zeroth log-moment, 
p' , only contributed to non hadronic LEC's correspond- 
ing to vacuum polarization diagrams, like the photon and 
graviton wave functions and the quadratic current quark 



mass corrections to the vacuum energy density. As such, 
this scale dependence triggered by dimensional transmu- 
tation in the SQM is innocuous since QED and quantum 
gravity also renormalize the photon and graviton wave 
functions respectively. 

If we assume a Spectral Quark Model a la Georgi- 
Manohar and consider a low energy expansion along the 
lines explained in the paper we would generate the scale 
dependent term 



C 



N r 1 



(i-g 2 A )(F^UF^ R U^) 



(4tt) 2 6 

[±-g] 

-{l-g 2 A f{{D^D v U) 



«(1 - g A )(F^UD"U^ + F*D^D"U) 
1 

V 

1 



+ -(l-g 2 A ) 2 ((D^D»Uf) 



(A.3) 



This terms arises from the photon wave function con- 
tribution and because of the additivity of the (1 — g A ) 
terms and the external spin one fields. This yields con- 
tributions to the LEC's involving pseudoscalars which 
depend on the arbitrary scale generated by dimensional 
transmutation. It is not clear whether this deficiency can 
be attributed to the g A 7^ 1 assumption or the spectral 
quark model construction. 
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